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We report charged-particle pair correlation analyses in the space of A(/> (azimuth) and Arj (pseudo- 
rapidity), for central Au + Au collisions at ^/sjvjv = 200 GeV in the STAR detector. The analysis 
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involves unlike-sign charge pairs and like-sign charge pairs, which are transformed into charge- 
dependent (CD) signals and charge-independent (CI) signals. We present detailed parameterizations 
of the data. A model featuring dense gluonic hot spots as first proposed by Van Hove predicts that 
the observables under investigation would have sensitivity to such a substructure should it occur, and 
the model also motivates selection of transverse momenta in the range 0.8 < pt < 2.0 GeV/c. Both 
CD and CI correlations of high statistical significance are observed, and possible interpretations are 
discussed. 

PACS numbers: 12.38Mh, 12.38Qk 



I. INTRODUCTION 

The Search for a Quark-Gluon Plasma (QGP) 
[l|, has been a high priority task at the Relativis- 
tic Heavy Ion Collider, RHIC H. Central Au + 
Au collisions at RHIC exceed [4j the initial energy 
density that is predicted by lattice Quantum Chro- 
modynamics (QCD) to be sufficient forproduction 
of QGP @. Van Hove and others [1,0,11 have pro- 
posed that bubbles localized in phase space (dense 
gluon-dominated hot spots) could be the sources of 
the final state hadrons from a QGP. Such structures 
would have smaller spatial dimensions than the re- 
gion of the fireball viewed in this mid rapidity ex- 
periment and the correlations resulting from these 
smaller structures might persist in the final state of 
the collision. The bubble hypothesis has motivated 
this study and the model described in Ref. @ has led 
to our selection of transverse momenta in the range 
0.8 < p t < 2.0 GeV/c. The Hanbury-Brown and 
Twiss (HBT) results demonstrate that for ^snn = 
200 GeV mid rapidity central Au + Au, when p t > 
0.8 GeV/c the average final state space geometry for 
pairs close in momentum is approximately describ- 
able by dimensions of around 2 fm Q . This should 
lead to observable modification to the At? A<p corre- 
lation. The present experimental analysis is model- 
independent and it probes correlations that could 
have a range of explanations. 

We present an analysis of charged particle pair 
correlations in two dimensions — A</> and Ar] - 
based on 2 million central Au + Au events observed 
in the STAR detector at ^JKw = 200 GeV [T3] 1 . 
The analysis leads to a multi-term correlation func- 
tion (Section III C and E) which fits the ArjAcf) dis- 
tribution well. It includes terms describing corre- 
lations known to be present: collective flow, reso- 
nance decays, and momentum and charge conserva- 
tion. Data cuts are applied to make track merg- 
ing effects, HBT correlations, and Coulomb effects 
negligible. Instrumental effects resulting from detec- 
tor characteristics are accounted for in the correla- 
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tion function. What remains are correlations whose 
origins are as yet unclear, and these are the main 
topic of this paper. We present high statistical pre- 
cision correlations which can provide a quantitative 
test of the bubble model || and other quantitative 
substructure models which may be developed. We 
also address possible jet phenomena. These preci- 
sion data could stimulate other new physics ideas as 
possible explanations of the observed correlations. 

This paper is organized as follows. Section II de- 
scribes the data utilized and its analysis. Section III 
describes finding a parameter set which fits the data 
well. Section IV presents and discusses the charge 
dependent (CD) and charge independent (CI) sig- 
nals. Section V presents and discusses net charge 
fluctuation suppression. Section VI discusses the 
systematic errors. Section VII is a discussion sec- 
tion. Section VIII contains the summary and con- 
clusions. Appendices contain explanatory details. 



II. DATA ANALYSIS 
A. Data Utilized 

The data reported here is the full sample of STAR 
events taken at RHIC during the 2001 running pe- 
riod for central Au + Au collisions at ^/snn = 200 
GeV. The data were taken using a central trigger 
with the full STAR magnetic field (0.5 Tesla). 

The central trigger requires a small signal, coinci- 
dent in time, in each of two Zero Degree Calorime- 
ters, which are positioned so as to intercept spec- 
tator neutrons, and also requires a large number of 
counters in the Central Trigger Barrel to fire. Ap- 
proximately 90% of the events are in the top 10% of 
the minimum bias multiplicity distribution, which is 
called the 0-10% centrality region. About 5% are 
in the 10-12% centrality region with the remainder 
mostly in the 12-15% centrality region. To investi- 
gate the sensitivity of our analyses to the centrality 
of the data sample fitted, we compared fits of the 
entire data sample to fits of the data in the 5-10% 
centrality region. Both sets of fits had consistent sig- 
nals within errors. Thus no significant sensitivity to 
the centrality was observed in the correlation data 
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from the central triggers. 

About half the data were taken with the mag- 
netic field parallel to the beam axis direction (z) 
and the other half in the reverse field direction in 
order to determine if directional biases are present. 
As discussed later in this subsection, our \ 2 analyses 
demonstrated there was no evidence of any difference 
in the data samples from the two field directions, and 
thus no evidence for directional biases. 

The track reconstruction for each field direction 
was done using the same reconstruction program. 
Tracks were required to have at least 23 hits in 
the TPC (which for STAR eliminates split tracks), 
and to have pseudo-rapidity, ij, between -1 and 1. 
Each event was required to have at least 100 pri- 
mary tracks. These are tracks that are consistent 
with the criteria that they are produced by a Au + 
Au beam-beam interaction. This cut rarely removed 
events. The surviving events totaled 833,000 for the 
forward field and 1.1 million for the reverse field. 
The transverse momentum selection 0.8 < pt < 2.0 
GcV/c was then applied. 

Based on the z (beam axis) position of the primary 
vertex the events were sorted into ten 5 cm wide bins 
covering -25 cm to +25 cm. The events for the same 
z bin, thus the same acceptance, were then merged 
to produce 20 files, one for each z bin for each sign 
of the magnetic field. 

The files were analyzed in two-dimensional (2-D) 
histograms of the difference in rj (Arj), and the dif- 
ference in 4> (A<f>) for all the track pairs in each event. 
Each 2-D histogram had 72 A<j) bins (5°) from -180° 
to 180° and 38 At] bins (0.1) from -1.9 to 1.9. The 
sign of the difference variable was chosen by labeling 
the positive charged track as the first of the pair for 
the unlikc-sign charge pairs, and the larger p t track 
as the first for the like-sign charge pairs. Our label- 
ing of the order of the tracks in a pair allows us to 
range over four A<p - Ar] quadrants, and investigate 
possible asymmetric systematic errors due to space, 
magnetic field direction, behavior of opposite charge 
tracks, and systematic errors dependent on p t . Our 
consistently satisfactory results for our extensive x 2 
tests of these quadrants for fits to this precision data 
revealed no evidence for such effects. 

Then we compared the A<p - Ar] data for the two 
field directions on a bin by bin basis. In the reverse 
field data, we reversed the track curvature due to the 
change in the field direction, and changed the sign of 
the z axis making the magnetic field be in the same 
direction as the positive z direction. This is done by 
reflecting along the z axis, and simultaneously re- 
flecting along the y axis. In the two dimensional A<fr 
- Ar] space this transformation is equivalent to a re- 
flection in A<j) and Ar). For each pair we changed the 
sign of its A(j> and A77 in the reverse field data. We 



then calculated a \ 2 based on the difference between 
the forward field and the reverse field, summing over 
the A<j) - Ar] histograms divided by the errors added 
in quadrature. The resulting x 2 f° r the t w0 fields 
showed an agreement to within 1.5a. Therefore we 
added the data for the two field directions. 

We also found there was no significant dependence 
within 2.2cr on the vertex z co-ordinate. Following 
the same methodology we added the files for those 
10 bins also. 



B. Analysis Method 

Separate A<p - Ar] histograms were made for 
unlikc-sign charge pairs and like-sign charge pairs 
from the same-event-pairs, since their characteristics 
were different. Both histograms are needed to later 
determine the CD and the CI correlations as defined 
in Sections IV A and B. These data are shown in Fig. 
1. Similar histograms were made with each track 
paired with tracks from a different event (mixed- 
event-pairs), adjacent in time, from the same z ver- 
tex bin. This allows the usual technique of divid- 
ing the histograms of the same-events-pairs by the 
histograms of the mixed-events-pairs which strongly 
suppresses instrumental effects such as acceptance 
etc., but leaves small residual effects. These include 
those due to time dependent efficiency of the track- 
ing in the readout boundary regions between the 12 
TPC sectors, which had small variations in space 
charge from event to event. Also space charge in 
general can cause track distortion and efficiency vari- 
ation event by event. The ratio of same-event-pairs 
to mixed-event-pairs histograms is shown in Fig. 2 
where the plot was normalized to a mean of 1. 

The expected symmetries in the data existed 
which allowed us to fold all four A</> - Ar] quadrants 
into the one quadrant where both A<j) and A77 were 
positive. After the cuts described in subsection C, 
we compared the unfolded bins to the folded aver- 
age for unlike-sign charge pairs and like-sign charge 
pairs separately. The \ 2 were less than the Degrees 
of Freedom (DOF), and the x 2 /DOF were within 
2-3cr of 1. 

Thus though we searched carefully in a number 
of ways to find asymmetries in the data via ex- 
tensive x 2 analyses and observation of fit behavior, 
none of any significance were found. By folding four 
quadrants into one, we quadrupled the statistics in 
each bin analyzed. Fig. 3 shows this folded data 
after cuts (see subsection C) were made to make 
track merging, Coulomb, and HBT effects negligi- 
ble. Henceforth the folded and cut data will be used 
for our fits. 
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Unlike-sign pairs 




Like-sign pairs 




FIG. 1: a) left side: The Acj> - Ar\ correlation data for unlike-sign charge particle pairs from the Star central trigger 
dataset shown in a 2-dimensional (2-D) perspective plot. The particle tracks have 0.8 GeV/c < pt < 2.0 GeV/c and 
r/| < 1.0. The structure that looks like tiles on a roof is due to the readout boundary effects of the 12 sector TPC. 
b) right side: The similar correlation data for like-sign charge particle pairs is shown. 



C. Cuts 

At small A<p and small Ai] (i.e. small space an- 
gles) track merging effects occur. To determine the 
cuts needed to reduce these effects to a negligible 
level, we varied small A(f> and small Ar] cuts. Simul- 
taneously the x 2 °f an approximate fit to the data 
using equation (3) + (4) + (5) was studied as a func- 
tion of the bins included in the fit. With larger cuts 
the x 2 behaved properly until one or more of the 
bins cut out was included in the fit. This caused a 
huge increase in x 2 , revealing that those bin(s) were 
distorted by merging effects. We confirmed by visual 
inspection that track merging effects clearly became 
important in the bins cut-out. The track merging 
caused a substantial reduction in track recognition 
efficiency, and supported the quantitative results of 
our x 2 analysis. The sharp sensitivity of the x 2 to 
the cuts showed how localized the merging losses 
were. We adopted these cuts and they are included 
in the folded data analyzed 

For the unlike-sign charge pairs we cut out the 
data from the regions 0.0 < Arj < 0.1 when 0° < 
A(f> < 20°, and 0.1 < Ary < 0.2 when 0° < A</>< 10°. 
These cuts to make track merging effects negligible 
also made the Coulomb effect negligible, since it is 
effective within a Ap of 30 MeV/c. This corresponds 



to a 2° opening angle at p t of 0.8 GeV/c. The 2° 
along the beam axis near 90° corresponds to a A-q 
of 0.035. 

For the like-sign charge pairs we cut out the data 
from the region 0.0 < Ai] < 0.2 when 0° < Acf> < 5°. 
The HBT effect is already small because of track 
merging within a Ap of 50 MeV/c, corresponding to 
a 4° opening angle at 0.8 GeV/c. The 4° along the 
beam corresponds to a A-q of 0.07. All the cuts de- 
scribed above were applied to the mixed events also. 
Thus the final track merging cuts selected reduce 
the tracking problems due to overlap and merging 
to negligible levels, and also reduce Coulomb and 
HBT effects to negligible levels. 

The fits to the subsequently shown data were 
made over the whole A(f> range. The data for 
|Ary| > 1.5 were cut out, since the statistics were 
low (see Fig. 1) and tracking efficiency varies at 
higher |A»7|. 



III. PARAMETERIZATION OF DATA. 

We want to obtain a set of functions which will 
fit the data well, and are interpretable, to the ex- 
tent practical. We utilized parameterizations rep- 
resenting known, expected physics, or attributable 
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FIG. 2: a) left side: The correlation data for the ratio of the histograms of same-event-pairs to mixed-event-pairs for 
unlike-sign charged pairs, shown in a two-dimensional (2-D) perspective plot A<f> - Arj. The plot was normalized to 
a mean of 1. 

b) right side: The similar correlation data for like-sign charge pairs. 



to instrumentation effects. Any remaining terms re- 
quired to obtain good fits to the data can be con- 
sidered as signals of new physical effects. Thus sig- 
nal = data - (known and expected) effects. The 
three known effects: elliptic flow, residual instru- 
mental effects, momentum and charge conservation 
terms were parameterized. We then find parameters 
for the signal terms which are necessary in order to 
achieve a good fit to our high statistical precision 
data. 

In the following we discuss the relative importance 
of the effects which are determined by x 2 m ou r final 
overall fit to equation (3) + (4) + (5). 



A. Elliptic Flow 

Elliptic flow is a significant contributor to two- 
particle correlations in heavy ion collisions, and must 
be accounted for in our analysis. Elliptic flow is not 
a significant contributor to two-particle correlations 
in p p collisions. The elliptic flow for STAR Au + 
Au data has been investigated in an extensive se- 
ries of measurements and analyses 11| using four 
particle cumulant methods and two particle correla- 
tions. This determines a range of amplitudes (v^) 
for cos(2A0) terms between the four particle cu- 



mulant which is the lower flow range boundary of 
W2(w2=0.035), and the average value of the four 
and two particle cumulant which is the upper range 
boundary of ^2(^2=0. 047). These results were deter- 
mined for the pt range 0.8 to 2 GeV/c by weighting 
the central trigger spectrum as a function of multi- 
plicity and centrality compared to the minimum bias 
trigger data for which the vi was measured. 

The constraint is used in the following way. We 
find our best fit to the data by allowing the parame- 
ters in the fit to vary freely until the \ 2 is minimized. 
However, there is in the case of the vi parameter a 
constraint that it must lie in the flow range cited 
above. Since our best fit was consistent with the 
bottom of the flow range it was not affected by the 
constraint. The value of vi cannot be varied to lower 
values because of the constraint and thus our best fit 
corresponds to the lower boundary of the flow range 
and had a vi of 0.035. However, vi can be varied to 
higher values until cither we reach the upper bound- 
ary or the fit worsens by lcr which is a \ 2 change of 
32. As discussed in Section VI, our fit worsens by 
lcr when v-i is increased to 0.042 well within the flow 
range 0.035 to 0.047. This represents a systematic 
error due to flow and is the dominant error( see Sec- 
tion VI ). In Appendix A, a detailed description of 
multi-parameter fitting and error range determina- 
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FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on A<f> and A77 based on demonstrated 
symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing 
with folded data only, b) right side: Similar folded data for like-sign charge pairs. 



tions are given. 



B. Instrumental Effects 

In addition to the tracking losses resulting from 
overlapping tracks in the TPC, which we handled 
by cuts, there are tracks lost in the boundary areas 
of the TPC between the 12 sectors. These regions 
result in a loss of acceptance in the <fi measurements 
since the particle tracks cannot be readout by elec- 
tronics. This produces a cf> dependence with a pe- 
riod of 30° which is greatly reduced from about a 
4% amplitude to about 0.02% (reduction by a factor 
of 200), by the normalization to mixed-event-pairs. 
However it is not completely eliminated, since STAR 
still has small time dependent variations from event 
to event, such as space charge effects and slight dif- 
ferences in detector and beam behavior. To correct 
for the residual readout A(f> dependence, a term with 
a 30° period to represent the TPC boundary peri- 
odicity and a first harmonic term with a 60° period 
were used. The unlike-sign and like-sign charge pairs 
behave differently over the TPC boundary regions in 
(j), because positive tracks are rotated in one direc- 
tion and negative tracks are rotated in the opposite 
direction. This requires that the sector terms have 
an independent phase associated with each. In the 



Tables I and II the terms are labeled "sector" and 
"sector2", "phase" and "phase2". The functional 
form for these sector effects is: 

Sectors(Z\</>) = S sm(12A(f> - tp) + S & sm{6A(j) - tp 2 ) 

(1) 

There is a <p independent effect which we attribute 
to losses in the larger r\ tracking in the TPC. We 
utilized mixed-event-pairs with a similar ^-vertex to 
take into account these losses. Imperfections in this 
procedure leave a small bump near Aij = 1.15 to be 
represented in the fit by the terms labeled etabump 
amp and etabump width in the tables. The width 
of this bump should be independent of the charge of 
the tracks, so we constrained it to be the same for 
like and unlike-sign charge pairs to improve the fit 
stability. Thus the functional form for this instru- 
mental effect is: 

Etabump(Z\r/) = Ee -(^-i-i5) 2 /2a% (2) 

In our final 3<7 fits for the combined like and 
unlike-sign charge pairs (See Tables I and II), leav- 
ing out the instrumental corrections would cause our 
3a fit to deteriorate to a 16a fit which is unaccept- 
able. It should be noted that standard x 2 analyses 
are often not considered credible if the fit exceeds 
2-3cr statistical significance. 
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C. Correlations associated with Momentum 
and Charge Conservation 

It is important to ensure that momentum and 
charge conservation correlation requirements are 
satisfied. For random emission of single particles 
with transverse momentum conservation globally 
imposed, a negative cos(A(f>) term alone can repre- 
sent this effect. However, the complex correlations 
that occur at RHIC result in multiple uncorrelated 
sources which are presently not understood. It was 
not possible to fit our data with the cos(A0) term 
alone. This fit was rejected by IOOct. This was 
not surprising since random emission of single parti- 
cles with momentum conservation would not lead to 
the particle correlations observed at RHIC. There- 
fore we suspected that a more complete description 
of the momentum and charge conservation was re- 
quired. No one has succeeded in solving this com- 
plex problem in closed form even in the theoretical 
case where you observe all particles. Hence a rea- 
sonable approach was to try to solve it for the tracks 
we are observing in order to obtain a good fit. We 
used Fourier expansion in the two variables we have, 
namely A(f> and Arj. 

Assuming that the cos(A0) term for random sin- 
gle particle emission was the first term in a Fourier 
expansion of odd terms, a second term cos(3A0) was 
added and found to account for about 95% of the 
lOOcr rejection. Based on the residual analysis we 
concluded the remaining 5% required Aij dependent 
terms for its removal. Therefore we multiplied terms 
of the type cos(A0) and cos(3A^) by a A?y depen- 
dent polynomial expansion cutting off at Ar/ 2 . This 
essentially removed the remaining 5% rejection. 

Some of the Ar/ dependent terms make quite small 
contributions in the fits, but in the aggregate they 
improve our overall fit by about 4cr. Thus a good 
fit of 3<T is downgraded to an unacceptable one of 
7cr without the A77 dependent expansion, since they 
add. One should note that fits which exceed 3cr are 
often not considered credible by experienced prac- 
titioners of precision data analysis, and a precision 
data analysis was a prime objection of this paper. 
In addition we found that we needed an overall A-q 2 
dependent fit parameter. 

If we take the sum of the terms we have in the 
above subsections A, B and C we obtain 



Name 



Value 



lump amplitude (A u ) 0.02426 , 



lump A(j) width (o"a</>) 



32.68T 



lump At; width (o"Atj) 



1.0581 



fourth (f) 



O.lOO; 



constant (Boo) 



0.99497; 



A V 2 (Bo 



0.00078=f0.00018 



cosA<^) (-B10) 



-0.00710; 



At; cos A0 (Bu) 



0.00092=f0.00078 



At? 2 cos A<f> (gig) -0.00049=f0.00057 



cos 3 A(f> (B 3 



0.00058 n 



0.00020 



At? cos 3 A0 (B31) 



0.00030±0.00062 



A77 2 cos 3A^ (B32) -0.00014=f0.00044 



sector (S) 



0.00016±0.00006 



phase (</?) 



,6±2.1 



sector2 (52 



0.00006±0.00007 



phase2 (^2) 



23.0±10.0 



A^bump amp (E) 0.00030=f0.00025 



AT/bump width (cte) 



0.189 (fixed) 



X 2 /DOF 



572 / 517 



TABLE I: Unlike-Sign Charge Pairs Fit Parameters for 
equation (3) + equation (4). The table has 3 sections. 
The top section lists parameter names, values and er- 
rors for the Approximate Gaussian Signal fit (lump). 
Fourth(f) is the additional term in the exponent. The 
source of the above is Subsection III E. The upper error 
is the change in each parameter when one increases the 
elliptic flow until the fit is la worse. The lower error is 
determined by varying each parameter at the lower range 
of the elliptic flow while all the other parameters are free 
to readjust until the fit loses la in significance. The mid- 
dle section has the normalization and a small background 
term followed by six Momentum and Charge Conserva- 
tion terms. Source is Subsection III C. The bottom sec- 
tion has instrumental terms. Four terms which are due to 
TPC sector gaps and two due to large r\ tracking errors. 
Source is Subsection III B. Uncertainties are dominantly 
systematic, assessed as described in Appendix A. 



(3) 



Bk = (Known + Expected)Effccts = 
■Boo + B 02 A77 2 + B 10 cos A(f> + Bu Arj cos A<p 
-B 12 Arj 2 cos A<p + 2v\ cos(2A^6) + B 30 cos(3A^>) 
+B 3 iA77cos(3A(?!)) + B 32 Atj 2 cos(3A0) 
+ Sectors ( A <f>) + Etabump(A?7) 



D. Fitting with Bk 

We used the well known result [l2| that for a 
large number of degrees of freedom (DOF), where 
the number of parameters is a small fraction of DOF 
and the statistics are high, the \ 2 distribution is nor- 
mally distributed about the DOF. The significance 
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Name 



Value 



lump amplitude (A t ) 0.01823 



-0.00482 
+0.00069 



lump A<p width (ua^O 



32.02 



-2.91 
-1.02 



lump A 7? width (cta^) 



1.847] 



dip amplitude (A d ) -0.00451 



+ 0.00090 



dip A(f) width (<7A0d) 



14.23 



-2.64 
+2.91 



dip Arj width (tJA-qd) 



0-2281°;°°° 



constant (Boo) 



0.995811°;°™! 



Ar? 2 (Boa) 



0.00100=F.00017 



cosA<^> (Bio) 



-0.007371 



1.00221 
'. 1)0040 



Ar;cos Acj) (Bn) 



0.00075 



+o 



00072 
00073 



A?7 2 cosA0 (B i2 ) 



-0.00015 



±11 



.00054 
00055 



cos3A0 (B30) 



0.000701° 



00033 
00020 



Ar? cos 3A0 (B 3 



-0.00027 



00064 
00065 



Arf cos3A0 (B32) 



0.00026 



00044 
00045 



sector (S) 



0.00021±0.00006 



phase (ip) 



22.6=Fl-5 



sector2 (S2) 



0.00007±0.00007 



phase2 (1P2) 



32.8=F7.0 



Ajybump amp (E) 



00022+°'°°° 24 

U.UUUZZ_ () 0Q022 



Arybump width (erg) 



0.189 (fixed) 



X / 



'DOF 



588 



519 



TABLE II: Like-Sign Charge Pairs Fit Parameters for 
equation (3) + equation (4). The table has 3 sections. 
The top section lists parameter names, values and errors 
for the Primary Gaussian Signal fit (lump). This is fol- 
lowed by a much smaller Gaussian (dip) . Source is Sub- 
section III E. The upper error is the change in each pa- 
rameter when one increases the elliptic flow until the fit is 
la worse. The lower error is determined by varying each 
parameter at the lower range of the elliptic flow while 
all the other parameters are free to readjust until the fit 
loses lcr in significance. The middle section has the nor- 
malization and a small background term, followed by six 
Momentum and Charge Conservation terms. Source is 
Subsection III C. The bottom section has instrumental 
terms. Four terms which are due to TPC sector gaps and 
two due to large 77 tracking errors. Source is Subsection 
III B. Uncertainties are dominantly systematic, assessed 
as described in Appendix A. 



of the fit dec reases by lcr whenever the x 2 increase 
is equal to y/2(DOF) which for our 517-519 DOF is 
equal to 32. Appendix A gives further details. 

If we fit the functional form of Bk (equation (3)) 
to both the unlike-sign charge pairs and the like-sign 
charge pairs, i.e. the whole data set, x 2 is about 
10,400 for about 1045 degrees of freedom. 

The standard deviation on 1045 DOF is about 46 
so the fit is rejected by around 200o\ We used many 
free parameters (15) in equation (3) yet additional 
unknown terms appear to be so sharply varying that 



it is clearly impossible for the Bk functional form 
to fit the data. A more detailed discussion of the 
statistical methods used in the data analyses, and 
the treatment of systematic and other errors in this 
paper is given in Appendix A. 



E. Signal Terms 

Many signal terms in physics are gaussian-like. 
We therefore tried fitting the signal data using two 
dimensional (2-D) gaussians or approximate gaus- 
sians. The unlike-sign charge pairs correlations were 
well fit (2a) , by adding to Bk an additional 2-D ap- 
proximate gaussian in Arj and A<f) (Fig. 4a) given 
by: 



Unlike — sign Signal = A u e 



(4) 

Considering the enormous improvement in fit 
quality afforded by the addition of this signal term, 
we conclude that this function in equation (4) pro- 
vides a compact analytic description of the signal 
component of the unlike-sign charged pairs corre- 
lation data. The fit was improved by the addi- 
tion of a term dependent on (A7/ 4 ) in the exponent 
(called "fourth" in Table I). The fit for the unlike- 
sign charge pairs signal data is shown in Fig. 4a. 
The unlike-sign charge pairs signal data correspond- 
ing to the fit (unlike-sign charge pairs data minus 
Bk ) is shown in Fig. 4b. 

The like-sign charge pairs data which also could 
not be fit by Bk alone, was well described (2a) when 
we added (see Fig. 3b) a positive 2-D gaussian and 
a small negative 2-D gaussian dip given by: 

Like - sign Signal = ^-(^V^W/Ki) 

+A de -( A * 2 / 2 ^+ A " 2 / 2 (5) 

Therefore, we conclude that equation (5) provides 
an efficient description of the signal component of 
the like-sign charge pairs data. The large signal re- 
ferred to in Table II as "lump" in the like-sign charge 
pairs correlation is a 2-D gaussian centered at the 
origin. It is accompanied by a small narrower 2- 
D gaussian "dip" (also centered at the origin) sub- 
tracted from it. The terms are labeled in the fits 
as: lump amplitude, dip amplitude, lump Ae6 width, 
lump Arj width, dip A(f> width, and dip Ar? width. 
Note that the volume of the dip is about 1.6% of the 
volume of the large lump signal. However, if we ne- 
glected to include this small dip, our fit deteriorates 
by 28a. This dip might be caused by suppression 
of like-sign charge particle pair emission from local- 
ized neutral sources such as gluons. The function fit 
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to the like-sign charged pairs signal data is shown 
in Fig. 5a. The like-sign charge pairs signal data 
(like-sign charge pairs data minus Bk) is shown in 
Fig. 5b. The uncertainties quoted throughout this 
paper corresponds to a change in \ 2 of Ax 2 = 32, 
rather than the more commonly used Ax 2 = 1 (see 
Appendix A for more details). 



F. Summary of Parameterizations 

Equation (3) + (4) yields a 2a fit for the unlikc- 
sign charge pairs correlation, and thus is an appro- 
priate and sufficient parameter set. 

Equation (3) + (5) yields a 2a fit for the like-sign 
charge pairs correlation, and thus is an appropriate 
and sufficient parameter set. 

Equation (3) + (4) + (5) yields a 3a fit for the en- 
tire data set, unlike-sign charge pairs correlation + 
like-sign charge pairs correlation, and thus is an ap- 
propriate and sufficient parameter set for the entire 
data set. 

The complex multi dimensional x 2 surface makes 
the \ 2 increase non linearly with the number of error 
ranges (la) given in the Tables I and II. Therefore, 
in order to determine the significance of a parameter 
or group of parameters one must fit without them, 
and determine by how many a the fit has worsened. 
Then one uses the normal distribution curve to de- 
termine the significance of the omitted parameter(s). 
The results of such an investigation were: 

1. Leave out cos(3A</>) and the fit worsens by 95a. 

2. Leave out dip and the fit worsens by 28a. 

3. Leave out Sector terms and the fit worsens by 
11.7a. 

4. Leave out 77 dependent cos(A</>) and cos(3A</>) 
terms and the fit worsens by 4.3cr. 

5. Leave out ?7bump terms and the fit worsens by 
1.54a. 

6. Leave out sector 2 and the fit worsens by lAla. 

Since our global fit with all the above terms is at 
the 3a level, we conclude that the above terms are 
necessary. The huge number of a we get for some of 
the parameters(s) are only to be taken as a qualita- 
tive indication of the need for the parameter(s) to fit 
our high precision data. Due to the uncertainties in 
determining the number of a far out on the tails of 
the normal distribution (i.e. > 10er), a quantitative 
interpretation of them cannot be made. 



IV. CD AND CI SIGNALS 

A. Charge Dependent (CD) Signal 

If we subtract the entire like-sign charge pairs cor- 
relation (equation (5)) from the unlike-sign charge 
pairs correlation (equation (4)) we obtain the CD 
correlation. However it is observed that the back- 
ground (Bk) of the two terms are close enough in 
value to cancel each other in the subtraction. Thus 
the CD signal is essentially the same as the entire 
CD correlation. In the extensively studied balance 
function [l3| for the correlation due to unlike-sign 
charge pairs which are emitted from the same space 
and time region, it is argued that the emission corre- 
lation of these pairs can be approximately estimated 
from the balance function. The balance function for 
these charged pairs is proportional to the unlike-sign 
charge pairs minus the like-sign charge pairs. There- 
fore, the CD signal which is approximately equal to 
the CD correlation is a qualitative measure of the 
emission correlation of unlike-sign charge pairs emit- 
ted from the same space and time region. In addi- 
tion to the approximations involved in the balance 
function, there is the modification to the CD signal 
discussed in Section IV D. 

Fig. 6a shows the fit to the CD signal. Fig. 6b 
shows the CD signal data that was fit. The signal 
form is a simple gaussian in both Ac/) and A77. The 
gaussian width in A</> is 28.3° ± 3.4° (0.49 ± 0.059 
rad) and the gaussian Aij width is 0.485 ± 0.054. 
Converting this pseudo-rapidity to a angle yields 
a width of 26.7° ± 3.0° (0.47 ± 0.052 rad). This cor- 
relation has the same angular range in A77 and Ac/). 
Therefore we observe pairs of opposite charged par- 
ticles emitted randomly in the r\ and direction with 
a correlation which has an average gaussian width of 
about 27.5° ± 3.2° (0.48 ± 0.056 rad) in the 9 angle 
corresponding to Ar) and also in the azimuth Ac/). 

B. Charge Independent (CI) Signal 

If we add the like-sign charge pairs signal to the 
unlike-sign charge pairs signal, we obtain the CI sig- 
nal. The CI signal fit shown in Fig. 7a displays the 
average structure of the correlated emitting sources. 
Fig. 7b shows the CI signal data which was fit by 
the analytic distribution shown in Fig. 7a. 

In order to get a good fit to the shape, both in 
A77 and Ac/>, we needed a more complicated form 
than one simple 2-D gaussian. The CI signal actually 
contains an approximate gaussian from the unlike- 
sign charge pairs signal plus two gaussians from the 
like-sign charge pairs signal. 

We want to obtain a measure of the effective Ac/) 
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Unlike-sign signal fit 



Signal Unlike-sign 



0.02 



0.01 




(a) %J 50 



0.02 



0.01 




FIG. 4: "(Color online)" a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from 
the fit to the normalized and folded unlike-sign charge pairs signal data. 

b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side. 



and A-q widths of the overall pattern of the CI signal. 
A good method for doing this is to compare the CI 
signal with a single gaussian which yields the same 
root mean square (RMS) values as the actual good 
fit described above involving two gaussians and an 
approximate gaussian. For Acj) that a is 32.0° ±0.6° 
(0.56 ± 0.01 rad). For A77 that a is 1.55 correspond- 
ing to an angle of 66. Ol^f (1.15 ± 0.02 rad). Thus 
the correlation is about twice as wide in A9, the an- 
gle corresponding to Afy, than in the angle corre- 
sponding to Acj) (see Fig. 7a). 



Another STAR measurement reports CD and CI 
correlations at ^/snn — 130 GeV/c [l4j]. 2 . The ma- 
jor differences between that previous analysis and 
the present work are the larger range in p t 0.15 < 
p t < 2.0 GeV/c and the lower statistical quality of 
the limited 130 GeV/c dataset. Although the anal- 
ysis in Ref. [l4| included low p t particles there is 
reasonable qualitative agreement with the present 
results. 



2 The CD correlation in Ref. [T3| is defined as like-sign pairs 
minus unlike-sign charge pairs which has opposite algebraic 
sign relative to the definition used in this paper 



C. Resonance Contribution 

To determine the maximum contribution and ef- 
fect resonances can have on the CI signal we con- 
sider the following. Resonances mainly decay into 
two particles. Therefore, neutral resonances which 
decay into charged pairs could be a partial source of 
the unlike-sign charge pairs in CD and CI signals. 
However, resonances do not add significant correla- 
tion to the like-sign charged pairs in either the CD 
or the CI signals. 

We expect that the final state particles from the 
resonance decay distribution will be symmetrical in 
the Acj> angle and the angle corresponding to Ar] 
since we do not expect any polarization mechanism 
that would disturb this expected symmetry. The 
fact that in the CD both the Acj) angle and the an- 
gle corresponding to Arj are approximately the same 
supports this. In this section we are evaluating the 
possible effects of resonance decay on the CI signal 
which definitely is observed to be elongated in the 
corresponding Ar/ angle by about a factor of 2 com- 
pared to the Acj) angle. The CI signal is always de- 
fined as: CI signal = unlike-sign charge pairs signal 
+ like-sign charge pairs signal. CD signal = unlike- 
sign charge pairs signal - like-sign charge pairs sig- 
nal. By subtracting the CD signal from the CI signal 
and rearranging terms in the equation one obtains 
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Like-sign signal fit 



Like-sign signal 



0.01 




(a) %J 50 




(b) 4 r ™ 



FIG. 5: "(Color online)" a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of 
equation (5). 

b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side. 



the following. CI signal 
charge pairs signal). 



CD signal + 2- (like-sign 



There is particular interest in the shape of the 
CI signal when comparisons are made to theoretical 
models e.g. HIJING and Ref. [|. The A<j) shape is 
about the same in the CI and CD signals. There- 
fore, we need only estimate the effect of resonances 
on the Arj width. By making the extremely unre- 
alistic assumption that the CD signal is composed 
entirely (100%) of resonances, we can estimate the 
maximum effect of resonances on the CI signal Arj 
width. We use the equivalent gaussian which has 
the same root mean square (RMS) widths for this 
calculation. From the data we determine that the 
Ar) width of the CI signal is increased by a 7% up- 
per limit. However if we use the result from the 
Appendix B that the resonance content of the CD 
is 20% or less, we must divide the 7% by 5 which 
results in an increase of Arj width of only about 2%. 
Either estimate of Arj width increase is inconsequen- 
tial, compared to the observed approximate factor of 
two elongation of the A77 width compared to the A(f> 
width. Details of these calculations are given in Sec- 
tion VI Systematic Errors and Appendix B. 



D. Modification of CI and CD Signals 

The CI and CD signals existing at the time of 
hadronization are changed by the continuing inter- 
action of the particles until kinetic freeze-out, when 
interactions cease. The interactions are expected to 
reduce the signals. Therefore we expect that the ob- 
served signals are less than those existing at the time 
of hadronization. 



V. NET CHARGE FLUCTUATION 
SUPPRESSION 

Net charge fluctuation suppression is an observed 
percentage reduction in the RMS width of the dis- 
tribution of the number of positive tracks minus the 
negative tracks plotted for each event, compared to 
the RMS width of a random distribution. 

If there are localized, uncharged bubbles of pre- 
dominantly gluons in a color singlet, when the bub- 
ble hadronizes the total charge coming from the bub- 
bles is very close to zero. Therefore if we are detect- 
ing an appreciable sample of such bubbles, we expect 
to see net charge fluctuation suppression. 

It should be noted that net charge fluctuation sup- 
pression can be deduced from the CD correlation 
given previously. Net charge fluctuation suppres- 
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CD (Difference) signal fit 



CD (Difference) signal 



0.01 




0.01 




FIG. 6: "(Color online)" a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate 
gaussian which is symmetrical. The average gaussian width is 27.5° ± 3.2° (0.48 ± 0.056 rad) both in A(j> and A77. 
Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs 
emitted randomly on 2-D r)4> directions. 

b) right side: The CD signal data corresponding to the adjacent signal fit on the left. 



sion previously analyzed at lower momenta has been 
found consistent with resonance decay [l5| . How- 
ever the present analysis has different characteris- 
tics. Furthermore Ref.0] has made specific estimates 
for net charge fluctuation suppression applicable to 
this experiment. We performed a charge difference 
analysis for tracks within cuts of 0.8 < pt < 2.0 
GeV/c and \rj\ < 0.75. This would allow a compari- 
son with Ref. [|. 

For each event we determined the difference of 
the positive tracks minus the negative tracks in our 
cuts. There was a net mean positive charge of 4.68 
± 0.009. The width of the charge difference distri- 
bution given by (RMS) was 11.149±0.017. To deter- 
mine the net charge fluctuation suppression we need 
to compare this width with the width of the appro- 
priate random distribution, which would have no net 
charge fluctuation suppression. However we must ar- 
range a slight bias toward a positive charge so that 
we end up with the same net mean positive charge 
as observed. If we now assign a random charge to 
each track with a slight bias toward being positive 
such that the mean net charge is also 4.68 ± 0.009, 
the width (RMS) becomes 11.865 ± 0.017. The per- 
centage difference in the widths which measures the 
net charge fluctuation suppression is 6.0% ± 0.2%. 



VI. SYSTEMATIC ERRORS 

Systematic errors were minimized using cuts and 
corrections. The cuts (see Section II C) were large 
enough to make contributions from track merging, 
Coulomb, and HBT effects negligible. 

Systematic checks utilized x 2 analyses which veri- 
fied that the experimental results did not depend on 
the magnetic field direction, the vertex z coordinate, 
or the folding procedures (see Sections II A and B). 
By cutting the track 77 at 1.0 we keep the systematic 
errors of the track angles below about 1° [la ]. 

In Section IV C we referred to a simulation in 
Ref. [1], which estimated that the background res- 
onance contribution to the CD correlation cannot 
be more than 20% (see Appendix B). From the res- 
onance calculations discussed in Appendix B the 
background resonance contribution to the unlike- 
sign charge pairs correlation would be reduced to 
5%. The details that justify this are: Our anal- 
yses measured an average of 19 unlike-sign charge 
signal pairs per event from the CD correlation sig- 
nal. The total average number of unlike-sign charge 
signal pairs is 88. Therefore with the extreme as- 
sumption that all of the CD correlation signal is due 
to background resonances, 21.6% of the unlike-sign 
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FIG. 7: "(Color online)" a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits 
of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a A(f> width of about 
32° (0.52 rad), and Ar/ width of about 66° (1.15 rad) which is about double the A(f> width, 
b) right side: The CI signal data corresponding to the adjacent fit on the left. 



charge signal pairs are due to resonances. However 
as discussed in Appendix B the estimated contribu- 
tion of resonances to the CD correlation is 20% or 
less. Therefore 21.6%/5 = 4.3% is the resonance 
content of the unlike-sign charge signal pairs, which 
we rounded to 5%. 

We also concluded via an upper limit calcula- 
tion with the above extremely unrealistic assump- 
tion that the shape of the CI signal could only have 
the measured width in A-q increased by about 7%. 
This is inconsequential compared to the factor of two 
increase in Ar] width compared to A(f> width. For 
the simulation discussed above, in which resonance 
backgrounds contribute 20% to the CD correlation, 
the increase in the CI correlation Arj width would 
be limited to about 2%. 

The CI signal and the CD signal are more phys- 
ically significant than the unlike-sign and like-sign 
charge pairs signals, since they are physically in- 
terpretable. The CI signal fit displays the average 
structure of the correlated emitting sources (Section 
IV B). The analysis suggests that the CD signal 
(Section IV A) is a qualitative representation of the 
emission correlation of the unlike-sign charge pairs 
emitted from the same space and time region. El- 
liptic flow contributions to the like and unlike-sign 
charge pairs signals are approximately equal and 



therefore cancel in the CD signal. Thus the CD sig- 
nal is not affected by uncertainty in the elliptic flow 
amplitude. 

In Fig. 8a we show the CI signal data with the 
lower flow range value of mean elliptic flow ampli- 
tude (weighted over the p t range of the data) used 
in the fit, namely, v% — 0.035. This is our best fit, in 
which the value of V2 is consistent with the lower flow 
range. Thus Fig. 8a is the same as Fig. 7b. In Fig. 
8b we show the CI signal data with our maximum 
amount of elliptic flow allowed. This value of elliptic 
flow causes our fit to be la worse, and corresponds to 
a x 2 increase of 32. This results in a mean weighted 
t>2 = 0.042. This value of elliptic flow lies in the de- 
termined flow range, but is smaller than the upper 
limit of 0.047 [ill- 

The upper error ranges in Tables 
I and II are determined by the effect of this change in 
the elliptic flow (0.042> v 2 > 0.035), since it is the 
dominant error which results in these range values. 
The main change in Fig. 8b compared to Fig. 8a 
is that the peak amplitude is reduced by 25%, and 
the 2 dimensional area is reduced by 34%. However 
the most significant CI signal parameters in compar- 
ing to theoretical models are those that measure the 
shape, namely the ratio of RMS A(j> and Ar/ widths, 
of a gaussian equivalent to the fit. For our best fit 
Acj) width = 32.0° ± 0.6° (0.56 ± 0.01 rad), and A77 
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FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the 
best x 2 (the same as Figure 7b) plotted as a 2-D perspective plot on A0 vs. A77. 

b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis 
(x 2 was worse by la). 



width = 1.55 which is equivalent to 66.0° q 1 ^ (1.15 
± 0.02 rad). The ratio Ary equivalent angle/ Acf> = 
2.06. For the case of the maximum elliptic flow value 
the determined widths were A(f> = 30.1° (0.53 rad), 
and At? = 1.375 (61.6°or 1.08 rad). The ratio Arj 
j A(j) = 2.05. Thus the important shape ratio has 
changed about 1%. 

Let us now address the errors due to contamina- 
tion by including secondary particles arising from 
weak decays and the interaction of anti-protons and 
other particles in the beam pipe and material near 
the beam pipe. These secondary particle back- 
grounds have been estimated to be about 10-15% 
[17] . In this analysis we are concerned mainly with 
the angles of the secondary particles relative to the 
primaries that survive our high p t cut, not their iden- 
tity or exact momentum magnitude. Our correla- 
tions almost entirely depend on angular measure- 
ments of A<j) and At]. In the range 0.8 < pt < 2.0 
GeV/c we have considered the behavior of weakly 
decaying particles and other non-primary particles 
which could satisfy our distance of closest approach 
to the primary vertex and p t cuts. Because high p t 
secondary particles are focused in the same direc- 
tion as the primaries, only a fraction of these parti- 
cles have sufficient change in angle that would cause 
an appreciable error in our correlation. Hence our 



estimate is a systematic error of about 4% due to 
secondary particles. 

Below we summarize our extensive discussion of 
systematic errors in this section. Track merging er- 
rors, Coulomb, and HBT effects were made negligi- 
ble by cutting out effected bins at small space an- 
gles. Instrumental errors were corrected for in the 
parameterization. Elliptic flow error is our dominant 
systematic error, due to the uncertainty that exists 
in the elliptic flow analyses [ll[. It should be noted 
that the size of the upper error range of our pa- 
rameter errors in Table I and II are determined by 
the maximum value of the elliptic flow we allowed 
(^2=0. 042) which corresponds to a x 2 change of 32 
(la). The sizeable change to the CI signal shown in 
Fig. 8b compared to Fig. 8a is due to this maxi- 
mum allowed value of elliptic flow. As shown above 
fortunately the shape ratio of the CI signal which 
is important for theoretical comparison is changed 
only by about 1%. When more accurate elliptic flow 
results become available it is a simple matter to in- 
sert these in the fits and reduce the resultant error. 
The next smaller systematic error is due to possible 
background resonance errors. The smallest system- 
atic error is due to secondary contamination. As 
discussed above none of these errors affect our im- 
portant physical conclusions significantly. 
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VII. DISCUSSION 

Highly significant correlations are observed for 
unlike-sign charge pairs and like-sign charge pairs, 
and consequently for charge dependent (CD) and 
charge independent (CI) signals. The CD signal is 
well described (with a la fit) by a symmetrical 2-D 
gaussian with an RMS width of about 30° in both 
A<f> and Arj. Conservative simulations in Ref. Q in- 
dicate that the contribution to this CD signal from 
background resonance decay is less than 20%. Sim- 
ulations (see Appendix B) estimate that the only 
significant change due to the background resonances 
is an increase in the CD amplitude of 20% or less. 
The CI signal is more complex and is the sum of 
the unlike-sign charge pairs signal and the like-sign 
charge pairs signal fits. Therefore it contains an ap- 
proximate gaussian from the unlike-sign charge pairs 
signal, and a large positive gaussian plus a small neg- 
ative gaussian (dip) from the like-sign charge pairs 
signal. The dip contains only 1.6% of the like-sign 
charge signal volume. This small dip, observed for 
the first time, has high significance in the fit. This 
feature is consistent with what one would expect for 
suppression of like-sign charge pair emission from a 
localized neutral source such as gluons. 

A 2-D gaussian which yields the same RMS widths 
as our overall CI signal, has a a along the A</> direc- 
tion of 32.0° ± 0.6° (0.56 ± 0.01 rad). However the 
Ar\ a is 1.55 corresponding to an angle of approxi- 
mately 66° (1.15 rad) and thus is twice as wide. 

The mean charge difference for tracks within the 
chosen analysis cuts is 4.68 ± 0.009 net positive 
charges, and the RMS variation of this quantity from 
event to event is 11.149 ± 0.017. A random charge 
assignment constrained to produce the same mean 
net charge has a larger width of 11.865 ± 0.017. 
The difference, 6.0%±0.2%, measures the net charge 
fluctuation suppression. 

The HIJING model produces jets in our p t range 
which are nearly symmetrical in A<f> and the angle 
corresponding to Ar] [l8|. The proper way to com- 
pare our data with HIJING is to compare our CI 
signal with the above mentioned HIJING CI corre- 
lation from jets. As shown in IV B, our CI signal is 
highly asymmetric since the angle corresponding to 
Arj is about twice the angle A<fi thus strongly con- 
tradicting these basic characteristics of HIJING jets 
correlations. 

Ref. [191] presents a detailed comparison of the 
STAR data presented here with a model based on a 
ring of localized bubbles emitting charged particles 
from the central fireball surface at kinetic freeze out. 
Good consistency with the STAR data presented in 
the present paper is report ed[H|. 



VIII. SUMMARY AND CONCLUSIONS 

We performed an experimental investigation of 
particle-pair correlations in A(f> and Arj using the 
main Time Projection Chamber of the STAR detec- 
tor at RHIC. We investigated central Au + Au col- 
lisions at ^Jsmn = 200 GeV, selecting tracks having 
transverse momenta 0.8 < p t < 2.0 GeV/c, and the 
central pseudo-rapidity region \rj\ < 1.0. The data 
sample consists of 2 million events, and the symme- 
tries of the data in A77 and A(f> allow four quadrants 
to be folded into one. The entire data set (unlike- 
sign charge pairs and like-sign charge pairs) was fit 
by a reasonably interpretable set of parameters, 17 
for the unlike-sign charge pairs and 19 for the like- 
sign charge pairs. These parameters are small in 
number compared to the total number of degrees 
of freedom, which was over 500 for each of the two 
types of pairs. Every fit reported here using these 
parameters was a good fit of 2 to 3 a or less. 

Section VI discusses systematic errors. From our 
analysis of the systematic errors we conclude that 
the important features of our data and the conclu- 
sions drawn have not been significantly affected by 
the systematic errors. 

Section VII discusses the results and model fits. 

This paper serves as an excellent vehicle for mak- 
ing detailed comparisons with, and testing of various 
theoretical models such as the Bubble Model Q and 
other relevant models. 
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APPENDIX A: MULTI-PARAMETER 
FITTING IN THE LARGE DOF REGION 
AND SYSTEMATIC UNCERTAINTIES 

Let us now consider the detail of multi-parameter 
data analysis in the large DOF region in terms of 
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change of x 2 for our 517-519 n for unlike-sign and 
like-sign charge pairs respectively. A lcr change in 
the significance of the individual fits require the 
change in x 2 of 32. The reader is referred to Ref.fl2| 
from which we quote: "For large n(DOF), the 
X 2 p. d.f. (probability density function) approaches a 
gaussian with mean = n and variance(cr) squared = 
2n." For n>50-100 this result has been considered 
applicable, and it remains applicable and becomes 
more accurate as n increases toward infinity. Thus 
for our 517-519 n, lcr = V2n =32. 

The statistical significance of any fit in this paper 
can be obtained by the following procedure: The 
number of cr's of fit = (x 2 - n)/32. The number of 
cr's refer to the normal distribution curve. 

For large DOF(bins - parameters) fluctuations oc- 
cur because of the many bins. When one fits the 
parameters, they will try and describe some of these 
fluctuations. Therefore we need to check whether 
the fluctuations in the data sample are large enough 
to significantly distort the parameter values. This 
has traditionally been done by using the confidence 
level tables vs. x 2 which allows a reasonable deter- 
mination of the fluctuations due to binning. The 
approximation (described above) we have used is an 
accurate extension of the confidence level tables. 

Let us consider our method of assigning system- 
atic error ranges to the parameters. Our objective 
is to obtain error ranges which are not likely to be 
exceeded if future independent data samples taken 
under similar conditions are obtained by repetition 
of the experiment by STAR or others. We want to 
avoid the confusion and uncertainty of apparent sig- 
nificant differences in parameters when the fit to new 
data does not differ significantly from a previous fit. 
We allow each parameter, one at a time, to be var- 
ied (increased and then decreased) in both directions 
while all the other parameters are free to readjust 
until the overall fit x 2 degrades in significance by 
lcr. This corresponds to an increase of x 2 of 32 for 
the unlike-sign and like-sign charged pairs fits. The 
X 2 surface has been observed, and x 2 increases very 
non linearly with small increases of the parameter 
beyond the error range. 

When one assumes that the parameter error is 
given by a x 2 change of 1 in the best fit [2(J, this 
is correct for the case where you know the underly- 
ing physics. The purpose of our use of parameters is 
not to determine their true values, since we do not 
know the underlying physics. We use the parame- 
ters to determine an analytic description of the data 
that fits it within 2-3a\ 

Let us address how one determines the significance 
of a parameter or a particular group of parameters 
in the case of this multi-parameter analysis. It is en- 
tirely incorrect to compare the error ranges in Table 



I and II with the parameter value. The only correct 
way to determine the significance of a parameter or 
group of parameters in the multi-parameter fit is to 
leave the parameter (s) out of the fit and refit with- 
out the parameter (s). One then determines by how 
many cr's the overall fit has degraded compared to a 
change of 32 for lcr. 

One should note that the upper error range on a 
parameter in Tables I and II represents the change 
in value of the parameter as we increase v% , our dom- 
inant error, from our best fit until the fit worsens by 
lcr (x 2 increases by 32). The other parameters are 
free to readjust during the foregoing variation. The 
lower error range is the result of changing the pa- 
rameter in the opposite direction of it in the upper 
range while the other parameters are free to read- 
just until the fit decreases in significance by lcr. The 
largest error which dominates the upper systematic 
error range determination is the variation of the vi 
error. However the constraint of V2 to the flow range, 
and the fact that our best fit is at the lower end of 
the range means that flow is not varied when we de- 
termine the lower error. Therefore the elimination 
of the large flow variation leads to generally smaller 
lower error ranges and thus observed asymmetry be- 
tween the upper and lower error ranges. 

Finally if one fixes any one of the parameters at 
the end of the error ranges with the correct vi as- 
signed (0.042 upper and 0.035 lower) and lets all the 
other parameters readjust, one will achieve an over- 
all fit which is only worse by lcr (a x 2 increase of 
32). 



APPENDIX B: RESONANCE 
CALCULATIONS 

In this section we discuss background resonances 
which have the following meaning. Background res- 
onances are resonances that only contribute to the 
unlike-sign charge pair correlation. Any resonance 
that is emitted as part of a correlated emitter will 
also have a like-sign charge correlation. This is be- 
cause any of the charges will have another particle 
of the same charge to correlate with which is also 
emitted by the correlated emitter. The background 
resonances we are discussing are not part of the cor- 
related emitters we want to study, but will add to 
the CD correlation. 

In Ref. [|[ a table of resonances and particles was 
considered in a thermal resonance gas model. The 
resonance content was adjusted to cause the 20% 
net charge fluctuation suppression observed at lower 
momentum in a STAR experiment for central Au + 
Au collisions. The result of this model calculation is 
a 20% contribution of resonances to the CD corre- 
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lation. The p t range of the analysis 0.8 < p t < 2.0 
GeV/c tends to have the symmetric pairs of par- 
ticles from resonance production decay into angles 
within the range of our cuts. However asymmetric 
decays only contribute one particle of the pair and 
thus do not add to the correlation. Therefore effec- 
tively only a part of the estimated resonance con- 
tribution affects the correlation measurement. Our 
simulations reveal that the only significant change 
that the background resonance contribution would 
introduce in the CD correlation is an increase in am- 



plitude of 20% or less. The CI signal = CD signal + 
2- (like-sign charge pairs signal). The resonance con- 
tribution to the like-sign charge pairs signal is small. 
The affect of resonances on the CI signal elongation 
in the Ar/ width is estimated by the model simula- 
tion to be a factor five lower than the extreme max- 
imum case previously assumed, that the entire CD 
correlation is composed of background resonances. 
Thus this would result in a At] width elongation of 
the CI signal of approximately 2% instead of the 
upper limit of 7%. 
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